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摘 要 和 轴 流 式 冷却 循环 泵 是 舰 船 动力 系统 中 的 关键 设备 之 一 , 泵 内 非 稳 态 流动 诱发 的 压力 脉动 是 影响 泵 振动 水 平 的 主 
要 水 力 因 素 . 运用 特定 设计 方法 对 某 比 转速 为 471 的 循环 泵 进行 了 设计 并 通过 实验 测量 进行 了 性 能 验证 . 基于 非 定常 数 
值 计 算 , 获得 了 不 同 工 况 下 循环 泵 多 个 截面 上 压力 脉动 系数 的 变化 规律 . 研究 表明 ; 叶轮 进口 压力 脉动 系数 沿 径 向 呈 递 增 
趋势 ,， 轮 缘 处 流体 激励 能 量 最 大 且 主 要 集中 于 叶 频 ; 叶轮 出 口 压力 脉动 系数 在 轮 缘 和 轮 载 处 取得 极 大 值 ， 主 要 激励 频率 为 
叶 频 及 其 高 次 谐 频 ; 导 叶 出 口 压 力 脉动 系数 沿 径 向 基本 一 致 , 主要 激励 频率 为 1/5 叶 频 ; 叶轮 出 口 诱发 的 叶 频 压力 脉动 是 
决定 泵 振动 水 平 的 主导 因素 . 研究 结论 为 循环 泵 减 振 技术 的 发 展 提供 了 一 定 的 理论 文 撑 . 
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Numerical Study on Pressure ee of Axial-fow Cooling 
u 


Abstract Axial-flow cooling ci ng pump is one of the key equipment in ship power system. 
The pressure pulsation ind jits unsteady flow is the main hydraulic factor affecting the pump 
vibration level. The oh with a specific speed of 471 was designed according to a special 
method and the perfo ce verification was carried out through theexperimental measurements. 
Based on unsteady numerical simulation, the variation of pressure pulsation coefficients in three 
sections of the pump under different operating conditions werdarialyzed. The results show that 
the pressure pulsation coefficient at the inlet of impeller i Cheases along the radial direction. The 


fluid excitation at the shroud is maximum and mainlyXdepending on the blade passing frequency 


( 户 PF)，The pressure pulsation coefficient at the IN and hub in impeller outlet section also has 
a maximum value and the main excitation fre wu y is fppr and its higher harmonic frequency. 
The pressure pulsation coefficient in guide vas utlet section is basically the same along the radial 
direction and the main excitation frequer MN 1/5 fppr. The greater pressure pulsation coefficient at 
fepris, the higher vibration will be. The research achievements are helpful to provide some reference 
for furthering the vibration reduction technology. 
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轴 流 式 冷却 循环 泵 (以 下 简称 循环 泵 ) 是 舰 船 动 ”的 需求 
力 系统 二 回路 中 的 关键 设备 之 一 ,主要 用 于 主 冷凝 噪声 辐射 水 平 是 影响 舰 船 综合 战斗 性 能 的 关键 
器 中 工 质 的 冷却 以 及 其 它 设备 的 冷却 和 冲洗 . 当 舰 。 因素 , 降低 噪声 源 的 噪声 强度 、 控 制 振动 噪声 的 传递 
船 处 于 倒车 、 停 车 或 低速 正 车 状态 下 ,循环 泵 运行 ”途径 是 现在 和 未 来 舰 船 噪声 控制 技术 研究 的 重点 和 
于 主动 旋转 工 况 以 满足 舰 船 动力 系统 安全 稳定 运行 。 发 展 方向 种. 循环 泵 主动 旋转 工 况 下 泵 内 非 稳 态 流 
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动 诱发 的 压力 脉动 是 影响 泵 体 振动 噪声 水 平 的 关键 
水 力 因素 ， 其 诱发 的 低频 段 噪声 信号 是 舰 船 的 主要 
目标 特征 信号 之 一 . 近年 来 , 我 国 和 欧美 等 发 达 国家 
之 间 在 循环 泵 减 振 降 噪 技术 水 平 上 的 差距 正在 逐步 
缩小 ， 如 何 进一步 降低 水 力 因 素 诱 发 的 振动 噪声 水 
平 是 目前 研究 的 重点 和 难点 问题 ， 其 关键 之 一 就 在 
于 深入 对 各 种 非 稳 态 流动 诱发 的 激励 特性 认识 ， 因 
此 对 循环 泵 内 非 稳 态 流 激 特性 的 研究 具有 十 分 重要 
的 学 术 和 应 用 价值 . 

由 于 循环 泵 结构 的 特殊 性 ， 对 其 内 部 非 稳 态 流 
激 特性 的 试验 测量 较为 困难 且 周 期 较 长 . 随 着 CFD 
技术 的 日 益 成 熟 ， 国内 外 学 者 普遍 采用 基于 非 定常 
数值 模拟 的 手段 研究 泵 内 流 激 特性 2-4。 施 卫 东 
等 回 在 考虑 叶 顶 间 阶 的 情况 下 对 某 轴 流 模型 泵 在 
不 同 采样 频率 和 采样 时 间 下 的 全 流 场 进行 了 非 定常 
数值 计算 和 分 析 , 研究 指出 应 选取 较 高 的 采样 频率 ， 
同时 建议 采用 12 倍 旋 转 周 期 的 采样 时 间 ; 黎 净 军 
TO a a 
压力 脉动 的 影响 ， ea 
处 压力 脉动 的 主 频 ; Tsukamoto H 等 [研究 发 狐 
频 及 其 高 次 谐 频 是 泰 内 主要 激励 频率 ， 肥 网线 
力 的 旋转 是 由 叶轮 和 导 叶 间 的 动静 下 次 姑 

以 上 研究 表明 ， 采用 合 诈 的 驶 六 妆 他 


拉 注 柏 球 下 水 力 屋 证 潜 半 这 守信 抽 be 
常规 轴 流 泵 有 较 大 的 区 别 ， 因 此 其 内 部 非 稳 态 流 激 
特性 必然 有 别 于 现 有 研究 结论 ， 

本 文 以 某 舰 船用 循环 泵 为 研究 对 象 ， 运 用 特定 
设计 方法 进行 设计 并 基于 实验 测量 开展 泵 性 能 验证 . 
在 此 基础 上 , 通过 全 流 道 非 定常 数值 计算 , 对 不 后 
下 全 下 双 个 地 上 讨 力 及 动 系 才 的 灾 化 机 
行 分 析 ， 初 步 揭示 了 该 型 泵 内 轮 缘 固 壁 和 轮 访 鱼 系 
上 流体 激励 的 主要 频率 及 影响 泵 振动 水 平 的 主导 诱 
因 ， 为 该 型 循环 泵 的 安全 稳定 运行 和 优化 提供 了 理 
论 支撑 和 发 展 方向 ， 

1 循环 泵 物理 模型 及 性 能 验证 
1.1 循环 泵 物理 模型 

循环 泵 的 名 义 比 转速 为 471， 相 关 设 计 参 数 为 : 
流量 330 m3/h, 扬程 5.1 m, 转速 1450 r/min, 设计 
工 况 效率 大 于 80%, 汽 蚀 比 转速 大 于 1050. 

循环 泵 叶轮 和 导 叶 的 水 力 设计 方法 参照 舰 船用 
轴 流 泵 设计 方法 开展 ， 其 设计 要 点 为 :叶轮 通 径 为 
200 mm, 叶轮 外 缘 的 贺 周 速度 小 于 等 于 15 m/s, 叶 
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轮 环 量 分 布 系数 沿 外 法 线 方向 先 增 加 后 减 小 ， 叶轮 
和 导 叶 的 叶片 数 分 别 为 5 和 11. 

循环 泵 为 卧 式 ,采用 水 平 进口 、90? 弯 管 出 口 型 
式 , 叶轮 旋转 中 心 线 至 泵 弯 管 出 口中 心 线 距离 为 300 
mm。 循环 泵 叶轮 和 导 叶 均 采 用 数控 加 工 方式 成 型 ， 
从 而 保证 了 数值 计算 对 象 和 实验 对 象 的 一 致 性 . 
1.2 循环 泵 性 能 验证 

为 验证 循环 泵 设计 方法 的 有 效 性 以 及 数值 计算 
方法 的 精确 性 ,构建 了 如 图 1 所 示 的 闭 式 回 路 实验 
系统 并 开展 了 相关 的 泵 性 能 测量 实验 ， 


图 1 实验 系统 
Fig. 1 Experimental setup 


系统 中 流量 的 测量 精度 为 0.3%， 泵 进出 口 静 压 
的 测量 精度 为 0.075%， 功 率 和 转速 的 测量 精度 为 
0.2%， 符 合 GB/ 工 3216-2005 国家 标准 中 1 级 精度 
要 求 . 能 量 性 能 焦 监 和 汽 蚀 性 能 实验 中 运行 工 况 的 
调节 分 别 基 求 系 统 出 口 阔 门 开 度 调 节 和 系统 真空 度 
调节 而 实现 。 

汗 是 泵 设计 工 况 下 实测 的 性 能 参数 为 ， 扬 程 


效率 80.1%， 汽 蚀 比 转速 1100， 满 足 设计 


能 指标 . 


2 数值 计算 


2.1 计算 区 域 及 网 格 划 分 

以 循环 泵 实际 尺寸 参数 为 依据 进行 三 维 流 道 造 

， 结 合 循环 泵 内 流 特点 ， 将 计算 区 域 分 别 划 分 为 
ee 叶轮 区 和 导 叶 出 口 区 三 个 区 域 , 同时 在 模 
型 泵 进出 口 各 顺延 四 倍 管 径 长 度 的 直 管 段 以 提高 计 
算 的 稳定 性 , 计算 区 域 如 图 2 所 示 。 

鉴于 循环 泵 叶轮 和 导 叶 内 流 道 的 三 维 扭曲 特征 ， 
采用 ICEM 软件 中 对 不 规则 区 域 具有 强 适 应 性 的 非 
结构 化 网 格 对 上 述 区 域 进行 网 格 划 分 。 以 泵 扬程 和 
水 力 效率 随 网 格 数 增加 是 否 趋 于 稳定 为 依据 进行 网 
格 无 关 性 验证 ， 最 终 确定 的 各 计算 区 域 网 格 数 分 别 
为 : 进口 区 428501、 叶 轮 区 1956019、 导 叶 出 口 区 
946152。 
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进口 区 


叶轮 区 ” 导 叶 出 口 区 


出 口 
图 2 计算 区 域 


Fig. 2 Three-dimensional computational domain 


2.2 数值 计算 方法 

计算 模型 进口 采用 速度 进口 边界 条 件 ， 来 流速 
度 由 计算 工 帝 点 流量 确定 ,出 口 采用 压力 出 口 边界 
条 件 。 依据 循 环 泵 过 流 部 件 实测 粗糙 度 设置 各 固 壁 


的 表面 粗糙 度 ， 固 壁 采 用 无 滑 移 条 件 ， 近 壁 区 采 有 


标准 壁面 函数 法 .采用 SIMPLEC 算法 实现 压 岂 
速度 之 间 的 耦合 ， 对 压力 项 采用 Standard 格式 吉 
离散 ， 其 它 各 项 均 采 用 一 阶 迎风 格式 进 称 稻 向， 
差 收敛 精度 设置 为 0.00001. 


定常 计算 中 采用 Realizab ' 锻 流 模型 封闭 
控制 方程 组 , 非 定 常 计算 中 条 于 Realizable k-e 


模型 的 DES 滑 流 模型 封闭 控制 方程 组 ， 其中, 采用 
Realizable k-s 模式 模拟 接近 固 壁 边界 的 区 域 和 涡流 
长 度 尺度 小 于 最 大 格 点 尺寸 的 区 域 ， 对 其 它 区 域 使 
用 LES 模式 进行 求解 加. 为 满足 Courant-Friedrich- 
Levy 稳定 准则 ， 时 间 步 长 设置 为 叶轮 旋转 1? 所 对 
应 的 时 间 ， 即 114.9 hs, 采样 时 间 取 为 12 倍 旋转 周 
期 所 对 应 的 时 间 。 XU 
2.3 非 定常 数值 计算 初始 值 验证 RS 

为 提高 计算 的 稳定 性 和 收敛 速度 ， 通 常 将 定常 
计算 结果 作为 后 续 非 定常 数值 计算 的 初始 值 ， 因 此 
准确 的 定常 数值 计算 结果 是 有 效 的 非 定常 数值 计算 
的 重要 前 提 . 在 设计 工 况 流量 (Qa) 下 , 定常 数值 计 
算 所 得 的 泵 扬程 为 5.28 m， 和 实验 值 之 间 的 误差 仅 
为 2.9%, 在 0.9 Qa 和 1.1 Qa 时 的 误差 均 在 3.5% 以 
内 ,表明 非 定常 数值 计算 初 场 的 精确 性 符合 要 求 . 
2.4 监测 点 设置 

为 综合 分 析 循环 泵 流 道内 主要 特征 位 置 处 的 压 
力 脉动 特性 ， 分 别 在 泵 叶轮 进口 、 出 口 和 导 时 出 口 
截面 上 设置 压力 脉动 监测 点 . 各 截面 上 从 轮 载 至 轮 
缘 沿 径 向 均 布 4 个 监测 点 ， 各 监测 点 的 编号 及 位 置 
如 图 3 所 示 . 
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图 3 监测 点 位 置 示意 图 


Fig. 3 Location of monitoring points 


3 计算 结果 分 析 
本 文中 采用 的 无 量 纲 参数 压力 脉动 系数 定义 


< 人 0 加 
27 


式 中 , p 为 瞬时 压力 值 ，5 为 采样 时 间 内 的 平均 压 
力 , p 为 液体 的 密度 , U 为 叶轮 外 缘 的 圆周 速度 。 
3.1 设计 工 况 下 压力 脉动 分 析 

图 4 为 设计 工 况 下 叶轮 进口 截面 上 各 监测 点 处 
压力 脉动 的 时 域 图 . 由 图 可 知 , 在 叶轮 进口 处 , 压力 
脉动 的 幅 值 沿 径 向 外 法 线 方 向 呈 递 增 趋 势 ， 因 此 叶 
和 
SA。 
“© 


0.000 0.025 0.050 0.075 0.100 0.125 0.150 0.175 0.200 
1t/s 


0.000 0.025 0.050 0.075 0.100 0.125 0.150 0.175 0.200 
t/s 
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0.000 0.025 0.050 0.075 0.100 0.125 0.150 0.175 0.200 
l/s 


0.000 0.025 0.050 0.075 0.100 0.125 0.150 0.175 2 
1/s 


图 4 1.0Qa A 
Fig. 4 Pressure signals in time domain of i RS 


points at 1.0Q@&q 


图 5 ep 
频 域 图 . 由 图 可 知 叶轮 着 昌 压力 脉动 的 特征 频 
率 为 1/5 叶 频 (24.2 Hz)、 叶 频 (121 Hz) 及 二 倍 叶 频 
(242 Hz)。 沿 径 向 外 法 线 方向 , 各 监测 点 1/5 叶 频 压 
力 脉 动 幅 值 基本 一 致 、 叶 频 压 力 脉动 幅 值 呈 快 速递 
增 趋势 、 二 倍 叶 频 压力 脉动 幅 值 呈 缓慢 递增 趋势 . 通 
过 比较 各 监测 点 不 同 特征 频率 处 压力 脉动 幅 值 的 大 
小 可 知 ， 叶 轮 进口 轮 缘 处 流体 激励 能 量 最 大 且 激励 、 
能 量 主要 集中 在 叶 频 ， 叶 轮 进 口 轮 慌 处 流体 激 

量 主要 集中 在 1/5 叶 频 . 流体 产生 的 激励 必 

轴 系 或 固 壁 向 外 界 传递 或 辐射 因此 叶轮 进口 处 影 
响 振 动 噪声 水 平 的 关键 水 力 因 素 是 叶轮 外 缘 处 叶 频 
压力 脉动 幅 值 的 大 小 、 

图 6 为 设计 工 况 下 叶轮 出 口 截面 上 各 监测 点 处 
压力 脉动 的 时 域 图 . 由 图 可 知 , 在 叶轮 出 口 处 , 压力 
脉动 的 幅 值 沿 径 向 外 法 线 方向 呈 先 减 小 后 增 大 的 变 
化 趋势 是 叶轮 轮 缘 和 轮 载 处 压力 脉动 幅 值 的 大 小 基 
本 相等 ， 因 此 叶轮 出 口 处 流体 对 轮 融 轴 系 和 轮 缘 固 
壁 的 激励 水 平 基 本 相当 . 

图 7 为 叶轮 出 口 截面 上 各 监测 点 处 压力 脉动 的 
频 域 图 。 由 图 可 知 ， 叶 轮 出 口 处 压力 脉动 的 特征 频 
率 为 1/5 叶 频 、 叶 频 及 其 高 次 谐 频 . 各 监测 点 1/5 叶 
频 处 压力 脉动 幅 值 的 变化 规律 和 叶轮 进口 一 致 ， 叶 


频 和 二 倍 叶 频 压力 脉动 幅 值 沿 径 向 外 法 线 方向 呈 先 
减 小 后 增加 的 变化 趋势 , 约 在 75% 半 径 (B3) 处 取得 
极 小 值 . 通过 比较 各 监测 点 不 同 特征 频率 压力 脉动 
幅 值 的 大 小 可 知 ， 叶 轮 出 口 处 流体 激励 能 量 主要 集 
中 于 叶 频 ， 其 次 为 二 倍 叶 频 ， 轮 缘 固 壁 处 叶 频 激励 
能 量 略 小 于 轮 载 处 ,而 其 叶 频 高 次 谐 频 激励 能 晤 则 
略 大 于 轮 载 处 .结合 压力 脉动 的 时 域 图 可 知 ， 叶 频 
及 其 高 次 谐 频 对 轮 缘 固 壁 和 轮 纹 轴 系 的 综合 激励 基 
本 相等 . 

由 上 述 分 析 可 知 ， 时 域 图 中 压力 脉动 的 幅 值 和 
频 域 图 中 主要 特征 频率 的 脉动 幅 值 沿 径 向 的 变化 规 
律 相 一 致 ， 因 此 循环 泵 导 叶 出 口 处 仅 给 出 如 图 8 所 


0.004 


0.003 


0.002 


Amplitude 


0.001 


0.000 
Al 


图 5 1.0Qa 工 况 


Fig. 5 Pressure S 
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图 6 1.08a 工 况 下 叶轮 出 口 各 监 到 


Fig. 6 Pressure signals in time domain of impell 


points at 1.0Qda 


示 的 压力 脉动 频 域 图 . 由 图 可 入 导 叶 的 整流 作 
用 下 ， 导 叶 出 口 处 压力 脉 过 显 减弱 且 沿 径 向 
基本 一 致 ， 其 主要 特征 频率 海 1/5 叶 频 和 2/5 叶 频 。 


通过 比较 各 监测 点 不 同 特征 频率 压力 脉动 幅 值 的 大 
小 可 知 ， 导 叶 出 口 处 非 稳 定 流动 诱发 泵 体 振动 的 主 
频 为 1/5 叶 频 且 对 轮 缘 固 壁 和 轮 载 轴 系 的 激励 水 平 
相等 


100 
200 300 


400 
已 
“eguency 


图 7 1.08a 工 况 下 叶轮 出 口 各 监测 点 处 压力 脉动 频 域 图 


Fig. 7 Pressure spectrum of impeller outlet points at 1.0QCda 


图 9 为 轮 载 上 各 监测 点 处 的 压力 脉动 频 域 图 ， 


结合 图 5、7、8 中 压力 脉动 幅 值 的 变化 趋势 可 知 , 循 
环 有 泵 流 道中 非 稳 定 流动 诱发 的 1/5 叶 频 激励 沿 轴 向 
基本 一 致 ， 叶 频 激励 主要 体现 在 叶轮 出 口 处 且 其 激 


励 能 量 远 大 于 其 它 特征 频率 产生 的 激励 ， 由 此 可 见 
泵 叶轮 出 口 处 诱发 的 叶 频 压力 脉动 是 诱发 泵 振动 的 
主导 水 力 因 素 。 
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图 8 1.08a 工 况 下 导 叶 出 口 各 监测 点 处 压力 脉动 频 域 图 
Fig. 8 Pressure Spectrum of guide-vane outlet points 
at 1.0Qd 
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图 9 1.08a 工 况 下 轮 载 各 监测 点 处 压力 脉动 频 域 图 


Fig. 9 Pressure spectrum of hub points at 1.0Qd 
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3.2 不 同 工 况 下 压力 脉动 对 比分 析 

图 10 为 不 同 工 况 下 各 监测 点 处 压力 脉动 系数 
的 均 方 根 误差 对 比 图 ， 图 中 :， 横 坐标 为 监测 点 编 
号 , 各 监测 点 对 应 的 三 个 柱状 图 从 左 至 右 分 别 对 应 
0.98aq、1.08q 和 1.1Qa 工 况 ; 纵 坐 标 为 压力 脉动 系 
数 的 均 方 根 误差 ,其 表达 式 如 下 : 


2_ (Cri -© 
RMSE=N\N= (2) 
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式 中 , Op 为 压力 脉动 系数 的 平均 值 , n 为 压力 脉动 
系数 的 计算 个 数 . RMSE 表征 瞬时 压力 脉动 系数 相 
对 于 平均 压力 脉动 系数 的 离散 水 平 ， 其 值 越 大 表示 
压力 脉动 越剧 烈 。 
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图 10 不 同 工 况 下 各 监测 点 处 压力 脉动 系数 的 均 方 根 误差 


and 1.1QCd 


Fig. 10 RMSE value of monitoring points at 0.9Qd， 


可 各， 卫 的 加 ,压力 和 
小 后 增加 且 在 设计 工 况 下 取得 极 小 明 循 环 泵 
设计 工 况 即 其 最 优 工 况 .综合 戏 掀 亦 辣 工 况 下 各 测 
点 处 压力 脉动 系数 均 方 根 误 大 小 可 知 ,循环 
泵 导 叶 出 口 处 压力 脉动 最 小 六 其次 为 叶轮 进口 处 , 叶 
轮 出 口 处 压力 脉动 最 剧烈 . 在 相同 工 况 下 , 各 轴 截 面 
上 压力 脉动 的 变化 趋势 和 前 面 的 论述 相 一 致 .在 叶 
轮 出 口 处 , 0.984 工 况 下 的 压力 脉动 显著 大 于 1.0Qa 
工 况 下 的 压力 脉动 ,其 原因 可 能 是 小 流量 工 况 下 泵 
内 剧烈 的 涡流 分 离 所 引起 的 ， 由 此 也 说 明 本 文 所 采 


用 的 设计 方法 存在 一 定 的 局 限 性 ， es 


度 来 看 ， 基于 该 设计 方法 所 得 的 循环 系 不 适 念 水 
于 小 流量 工 况 . 
4 认 

通过 对 特定 设计 方法 下 循环 泵 压力 脉动 特性 的 
数值 计算 和 分 析 , 得 到 以 下 结论 : 

1) 叶轮 进口 轮 缘 处 流体 激励 能 基 最 大 且 激 励 能 
基 主 要 集中 在 叶 频 ， 轮 载 处 流体 激励 能 量 主要 集中 
在 1/5 叶 频 ， 

2) 叶轮 出 口 处 叶 频 和 二 倍 叶 频 压力 脉动 幅 值 沿 
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径 向 外 法 线 方向 呈 先 减 小 后 增加 的 变化 趋势 ， 叶 频 
及 其 高 次 谐 频 对 轮 缘 固 壁 和 轮 纹 轴 系 的 综合 激励 基 
本 相等 ,叶轮 出 口 处 的 叶 频 压力 脉动 水 平 是 决定 泵 
振动 水 平 的 主因 ， 

3) 导 叶 出 口 处 压力 脉动 系数 沿 径 向 基本 一 致 ， 
主要 激励 频率 为 1/5 叶 频 . 

4) 基于 该 设计 方法 的 循环 泵 不 适宜 运行 于 小 流 
量 工 况 。 
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